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Recently we found that, among recursive digital filters using satu- 
ration arithmetic to contend with overflow, a fundamental difference 
exists between second and higher order filters: the latter may sustain 
large-amplitude overflow oscillations. In this paper we have derived 
a new criterion expressly designed for determining when a given 
high-order recursive system using saturation arithmetic is free of ov- 
erflow oscillations. The new criterion, which is easy to use, follows from 
this result: we associate with the given system two trigonometric 
polynomials in 6 of degree equal to the order of the given system; if any 
linear combination of the polynomials with nonnegative weights is 
positive for all 9 in [0,ir] then the system is free of all nontrivial periodic 
oscillations. We prove that the new criterion subsumes certain well- 
known criteria, such as Tsypkin's criterion, from the literature on 
nonlinear systems. To illustrate, three classes of special systems are 
investigated, and in each case the new criterion gives substantial im- 
provements. Finally, the new test is applied in the synthesis of high 
order sections for a realistic eighth-order system. 



I. INTRODUCTION 

Recently 1 we made the unexpected observation that, among recursive 
filters employing saturation arithmetic, a fundamental difference exists 
between second and higher order filters, namely, the latter may sustain 
large-amplitude overflow oscillations. This observation proved to be 
timely since it coincided with the awareness that economies of scale 
coupled with various recent developments make highly attractive the 
use of high-order sections in filter realizations. It has also come to light 
that the problem of possible large-scale oscillations is of interest not only 
in data filtering but in other areas where the natural structure is a 
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high-order recursive system, e.g., code converters (DPCM -*• PCM) and 
speech synthesizers. 

The economies of scale derive from the fact that the overflow detection 
and correction circuits, an expensive part of present-day filters, are as 
many as the number of sections employed; thus if a realization is com- 
posed of fourth-order sections rather than the conventional second-order 
sections, then the number of such circuits may be expected to be halved. 
The recent developments alluded to earlier refer to the almost simul- 
taneous developments of inexpensive, lower-power-consuming semi- 
conductor read-only memories, and the concept of distributed arithmetic 
blocks 2 - 3 in which ROMs are used to implement digital filters. In a pi- 
oneering study R. B. Kieburtz 4 recently estimated that in a particular 
application a saving of about 30 percent in parts may be achieved over 
the conventional design through the use of fourth-order sections using 
saturation arithmetic and implemented by ROMs. 

Thus there is much to be gained if high-order sections can be used, 
and for this to happen it is first necessary to ensure that the highly de- 
structive overflow oscillations are not present in a particular design. It 
is apparent that there is a useful role for an effective criterion for de- 
lineating stable systems employing saturation arithmetic. It is possible 
to conceive of the situation where such a criterion is incorporated in the 
early design, i.e., the criterion is introduced as a constraint in the ap- 
proximation problem. The other possibility of an arithmetic different 
from saturation arithmetic to contend with overflow is not pursued in 
this paper. 

There do exist many such criteria in the literature on the stability of 
a class of nonlinear feedback systems (i.e., the systems in the Lurie 
problem 24 ) of which the one under consideration here is a member 5-11 ; 
the reader may consult Ref. 8 for a comparative evaluation of some of 
these criteria. These criteria are in some sense generalizations of Ny- 
quist's criterion for linear feedback systems. The reader will find in Sec. 
5.2 a statement, in the context of the present problem, of Tsypkin's 
criterion and the discrete circle criterion, two well-known examples of 
such criteria. Unfortunately it is known that when the nonlinearity in 
the system is the one associated with saturation arithmetic then these 
criteria, including Barkin's criterion, 7 - 8 are of limited utility since they 
are excessively pessimistic. Examples to this effect may be found here. 
Also telling is the fact that these criteria do not predict that all second- 
order systems using saturation arithmetic are free of oscillations, 8 a fact 
proven in Refs. 12-14 by arguments special to second-order systems. This 
is not totally unexpected in view of the fact that the systems-theoretical 
criteria apply to large classes of systems and nonlinearities and, con- 
comitantly, use relatively little information (restricted to the sector in- 
formation, symmetry, and monotonicity) about the nonlinearity. 
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Fig. 1 — Schematic of unforced high-order filter employing saturation arithmetic. 



In this paper we have derived new criteria expressly designed for the 
system employing saturation arithmetic. Underlying the new criteria 
is the observation that certain unique passivity conditions are operative 
in the case of saturation arithmetic. Both the observations regarding the 
passivity conditions as well as the technique we use for deriving the 
criterion are believed to be new. The main ingredient in the derivation 
is the observation that the expressions associated with the passivity 
conditions in periodic solutions possess remarkable structure; namely, 
they are quadratic forms involving circulant matrices. 

The systems considered in this paper are of the form (see Figs. 1 and 
2) 



;„ = F( L ajx n -j\ n = 0,1,. 



(1) 




Fig. 2 — The saturation arithmetic nonlinearity. 
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where m is the order of the system, (a/) are the coefficients, and F(«) is 
a nonlinear function associated with saturation arithmetic, namely 

F{a) = a if|<r|^l 

= sgn o-if |<r| > 1 (2) 

It is tacitly assumed that the underlying linear system in eq. (1) is ab- 
solutely stable, i.e., 

m 

\ m - Z aj\ m -J 9* for all \\\>l. (3) 

Thus any nontrivial solution of eq. (1) will necessarily have either 1 or 
—1 as an element and consequently such solutions are referred to as 
overflow oscillations. Note that we are following convention in ignoring 
quantization effects in the description of the filter in eq. (1); in investi- 
gations of large-scale oscillations it is natural to focus on the gross non- 
linearity. 

The main result of this paper (Sec. 2.4) is simply stated: for a given 
system of order m with coefficients {a,-}, we associate two polynomials 
of degree m in cos 0, namely, 



pi(6) = 1 - cos - £ o,|cos j6 - cos 0' -1)0} (4) 



and 



p 2 (0) = 1 + cos - E ajioMJB + cos 0" - 1)0}- (5) 

/-i 

If any linear combination of the polynomials with nonnegative weights 
is positive for all in [0,tt] then the system in eq. (1) with arbitrary initial 
conditions does not admit any nontrivial periodic solutions. Certain 
generalizations of this criterion are derived in Sec. VI. 

In Sec. 2.5 ("How To Use The New Test") we show that the criterion 
may be used in a straightforward manner by one of two methods. The 
first method calls for plots of pi(0), p 2 (0) and Pi(0)/p2(0) for in [0,tt]. 
The second requires the consistency of a set of linear inequalities to be 
checked. The second method may also be used for the generalized cri- 
terion in Sec. VI. 

In Sec. Ill we examine three classes of special systems in detail. The 
results for the following canonical example 1 are typical: in a fourth-order 
system the poles are taken to be all real and repeated at p where | p \ < 
1. For \p\ > 0.669 overflow oscillations are proven to exist. Tsypkin's 
criterion and the circle criterion guarantee the absence of oscillations 
for \p\ < 0.384. The new test guarantees the absence of oscillations for 
\p\ < 0.610. 

1682 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1977 



In Sec. IV we apply the test to an eighth-order filter in the TDM-FDM 
translator, an extensively studied application of digital filtering. We find 
that this system can be lumped into two fourth-order sections employing 
saturation arithmetic, neither of which can sustain overflow oscillations. 
Both sections fail Tsypkin's criterion and the circle criterion tests. 

Finally, in Sec. V we prove that the new criterion (i) easily gives the 
well-known result that overflow oscillations do not exist in second-order 
sections, and (ii) subsumes Tsypkin's criterion. 



II. THE CRITERION 
2.1 Passivity properties 

On account of the special form of the nonlinearity F the system in eq. 
(1) possesses certain simply stated but important properties which we 
interpret as passivity properties. The criterion we derive is a direct 
consequence of these properties. 

We may write eq. (1) as a linear system with a forcing sequence present 
by defining 

•» a f(£*J**-j) "(l¥»-i). n = 0,l,... (6) 

so that, from eq. (1), 

m 

x n - E OjXn-j + en, n- 0,1, ... (7) 

/'-I 

The above nonhomogeneous linear recursion is used throughout the 

paper. Our procedure is to translate the important features of F(-) into 

tractable constraints on \e n j. 

Every solution of eq. (1) possesses the following properties. 

Proposition 1 : 

e„(x n -x„-i)<0, n-lA... (8) 

and 

e n (x n + x n -i)<0, n = l,2,... (9) 

For proof we observe that if 1 2 ajx„-j | < 1 then e n = and both condi- 
tions are obviously valid. If 

m 

£ a/Xn-j > 1 then x n = 1 from (1) 

7=1 

and e n < from (6) 
and, of course, \x n -i\ < 1. Thus, in this case eqs. (8) and (9) are valid. 
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Similarly, if 2)a/x ra _; < — 1 then x n = — 1, e n > and as |x„_i| < 1, the 
relations remain valid. 

Equation (8) may be viewed as stating that the forcing term e n has 
opposite sign from (x n — x n -{) which we may interpret as the discrete 
analog of "velocity." Thus locally the forcing term acts to reduce the 
velocity. Similarly interpreting (x n + x n _i)/2 as the local "distance," 
eq. (9) states that locally the forcing term also acts to reduce the distance. 
For these reasons we view eqs. (8) and (9) as passivity properties. 

Note that the two above conditions imply the following weaker con- 
dition: 

e n x n < 0, n = 1,2, . . . (10) 

We find upon reflection that the latter condition is completely equivalent 
to the nonlinearity F in eq. (1) lying in the sector bounded by lines of 
slopes and 1: 

< F{a)/a < 1 for all a. (11) 

This sector information is exploited in various criteria 5-11 but the ad- 
ditional information in eqs. (8) and (9) is not. 

2.2 Equations for an oscillatory solution 

We state the equations associated with every oscillatory solution of 
period N of the nonhomogeneous recursion in eq. (7). We find that the 
equations in matrix form involve a circulant matrix. We put into per- 
spective some well-known results on circulants which are assembled in 
Appendix A. 

A periodic solution of eq (7) with period N has associated with it the 
following set of N equations involving the coefficients {a,}, the elements 
of the solution Xi,X 2 , • • • ,Xn and the corresponding forcing terms 
Ei,E 2 , • . . ,En- 

X x = cliXn + a 2 Xfli-i + h a m XN- m +i + ^i 

X 2 = fli-^i + a 2 X^ H h a m X/v-m+2 + E 2 



(12) 



Xn - o.iXn-1 + 02-Xn-2 H \- a m X^- m + E^ 

In matrix form, 

MX = E. (13) 

The interesting feature of the N XN matrix M is that it is a circulant 15-17 
since it is the following polynomial in the primitive N X N circulant 
matrix P (see Appendix A for definition of P): 

M = I - f aj-pJ (14) 
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The circulant matrices have been extensively studied in the past and we 
are in the fortunate position of knowing a great deal of their eigens- 
tructure.t In particular, the eigenvalues of M are 

m 

1 - L aje-iJ^/N k = 1,2, ...JSI (15) 

The eigenvectors of circulants are also known. The following remarkable 
property of the eigenvectors of circulants is of utmost importance in the 
paper (see next section): all N X N circulants have an identical set of 
eigenvectors, i.e., the eigenvectors do not depend on the constituents of 
the matrix. Thus, the eigenvectors of any N X N circulant are 1 : 

Uk = -J— i e -ik(2ic/N) > e -i2k(2w/N) t e -i(N-l)k(2ir/N) t lj* (lg) 

k = 1,2, . . . ,N. Thus the real and imaginary components of the elements 
of each eigenvector are sequences of equiseparated samples of a sine 
function. Although these are complex vectors and circulants are not 
generally symmetric, the eigenvectors of circulants share an important 
property with eigenvectors of symmetric matrices in that they form an 
orthonormal set, i.e., 

ujuj = hi- (17) 

In matrix notation, 

U*U = I (18) 

where the eigenvectors \uk ) have been arranged as columns of the matrix 
U. 



2.3 Another representation of the passivity properties 

We combine the above information with the passivity properties stated 
in Proposition 1 to obtain a compact and useful representation of the 
passivity properties that are valid if an oscillatory solution to (1) exists. 
As in the preceding section an oscillatory solution is assumed to be of 
period N. 

Note that we may write 

(X N ,X lt . . . ,X N -t) = {X l ,X 2 , . . . ,X N )P' = X'P' (19) 



f Recently we have had another occasion 18 to use the eigenstructure of the matrix M. 
Willson 19 investigates the matrix M from a different angle. 

1 We denote the conjugate transpose by the superscript *. In the case of real matrices it 
is also denoted by the superscript '. 
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where P is the primitive N X N circulant. Thus 

ZE n (X n -X n - 1 ) = X'(I-P')E 

n=l 



= X'(I - P')MX, from (13) 

(20) 



[m ml 

I - £ ajPJ - P"" 1 + £ ajPJ- 1 X 

where in the final step we have used, in addition to the expression for 
M, the relations P' = P^-i and P^ = I. The key observation about eq. 
(20) is that the matrix there, being a polynomial in P, is a circulant. 

We undertake a convenient change of coordinates to diagonalize the 
matrix in eq. (20). Let 

Z == U*X (21) 

where U is, as in Sec. 2.2, the unitary matrix of eigenvectors ofNxN 
circulants. Denoting the known eigenvalues (see Appendix A) of the 
matrix in eq. (20) by /%, k = 1,. . . ,N, we obtain 



£ E n (X n - X n _!) = £ \Z k 1 2 Re ixk (22) 



Now, 



Re ii k = 1 - cos \k(2ir/N)\ - £ a; [cos \jk2-K/N\ - cos \{j - l)k2Tr/N)\ 

7=1 

(23) 

To put eqs. (22) and (23) into the most convenient form, define the 
polynomial pi(-) where 

m 

pi(0) = 1 - cos 6 - £ aj{cosj6 - cos 0' - 1)0} (24) 

y=i 

We then have 

E E n (X n - X„-x) = E \Z k \* Pl (k2T/N). (25) 

n=l k=l 

We proceed in identical fashion to obtain a similar expression corre- 
sponding to the other passivity condition in Proposition 1. Note that 

N r m 

E E n (X n + Xn-r) = X' I - E ajPJ + P"-i 

n=l L j=l 

- E ajV- 1 ] X (26) 

7=1 J 

Because of the previously mentioned (and crucial) property that all N 
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X N circulants have identical sets of eigenvectors, the diagonalizing 
transformation is same as the one undertaken previously in eq. (21). 
Hence 

E E n (X n + X n - l )= E |Z*| 2 ReX, (27) 

n=l k=l 

where we have denoted the eigenvalues of the matrix in eq. (26) by (A*}. 
Here 

m 

Re \ k = 1 + cos (fe27r/iV| - £ a; [cos \jk2ir/N\ + cos {(/' - l)k2Tr/N\] 

7=1 

(28) 
k = 1,2, . . . , N. Thus for the final form we obtain 

E E n (X n + X n - X ) = L \Z k \2p 2 (k2ir/N) (29) 

n=l k=l 

where the polynomial p 2 (-) is defined to be 

m 

p 2 (6) = 1 + cos - E a ; {cos ;0 + cos 0' - 1)0) (30) 

Now certainly Proposition 1 implies that 

E E n (X n - X„-i) ^0 and E £*(** + X„-i) < (31) 

n=l n=l 

The above, together with eqs. (25) and (27), yields: 

Proposition 2: If a periodic solution of period N with elements 
(Xi,X 2 , . . . ,Xn) exists for the recursion in eq. (1) then 

E \Z k \^ Pl (k2ir/N)^0 (32) 

and 

E |Z fe |2p 2 (&27r//V)<0 (33) 

where Z, as given in eq. (21), is a transform of X and pi(0) and p 2 (0), given 
in eqs. (24) and (30), are two polynomials in cos of degree equal to the 
order of the system of eq. (1). 

For a fourth-order system (m = 4) the two polynomials are 

Pi(0) = (1 + ai) — (1 + ai — a 2 ) cos — (a 2 — a 3 ) cos 20 

— (as — 04) cos 30 — CI4 cos 40 (34) 
and 

p 2 (0) = (1 — cti) + (1 — ai — a 2 ) cos — (a 2 + 03) cos 20 

— (as + 04) cos 30 — 04 cos 40 (35) 
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The polynomials for second- and third-order systems are obtained from 
the above by setting 03 = 04 = and 04 = 0, respectively. 

In Fig. 3a and b we have plotted pi(0) and p 2 {6) for a particular 
fourth-order system. 

2.4 The main result 

It is only a short step from Proposition 2 to the main result which 
is 

Theorem 1: If for any «i ^ and <x 2 > 0, 

aiPi(0) + «2P 2 (0) > for all 6 in [0,x] (36) 

then nontrivial periodic oscillations do not exist as solutions to eq. (1). 

Proof: The proof is by contradiction. Suppose a nontrivial (i.e., X ^ 0) 
periodic solution of period N exists and also that the hypothesis of the 
theorem is valid. Then for such a solution 

«i L \Z k \2 pi (k2r/N) + a 2 £ \Z k \*p 2 (k2-K/N) 
k=\ k=i 

= £ \Z k \z\a lPl (k2T/N) + a2 p 2 {k2-x/N)\ 

> (37) 

from the hypothesis. However, from the passivity conditions summarized 
in Proposition 2, 

«i E \Z k 1 2 Pi(k2T/N) + « 2 E \Z k \ 2 p 2 (k27r/N) < (38) 
k=i k=\ 

which is a contradiction. QED. 

Note that if it is desirable to know only that oscillations of a particular 
period N do not exist for eq. (1) then the following is a sufficient condi- 
tion: 

There exist 

«i ^ 0, a 2 ^ such that a l p 1 (k2Tr/N) + ci2p 2 {k2ir/N) > (39) 
for k = 1,2, . . . ,N. 

2.5 How to use the new test 

Given an mth-order system, there are two simple and straightforward 
ways in which the above result may be used to determine whether the 
system does not admit overflow oscillations. 

The first method requires Pi(0), P2(0)» and Pi{0)/p 2 (d) to be plotted 
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Fig. 3— Plots of polynomials p\(0), p 2 (0), and otiPi(O) + ot2P2(0) for the following 
fourth-order system: a i = 1.1015710, a 2 = -1.6571120, a 3 = 0.7733805, a 4 = -0.45135546. 
In (c), a, = .92348761, a 2 = .16965636 

for 9 in the interval [0,tt] . When such plots are available the first step is 
to see if there is any 6 for which both pi(6) and p 2 (0) are negative; if this 
is the case the hypothesis of Theorem 1 obviously cannot be satisfied and 
the test is automatically failed. Assuming that this is not the case we find 
upon reflection that the hypothesis of Theorem 1 is satisfied if and only 
if 

max [pi(d)/p 2 (0)) < min [pi(0)/p 2 (0)] (40) 

|0|pi(0)>O,p 2 (0)<O| |0|pi(0)*SO,p2(0)>O| 

In fact, if the above is true the interval defined by the left- and right-hand 
sides of (40) is not empty and the hypothesis of Theorem 1 is satisfied 
by taking ai = 1 and -a 2 to be any value in the interval. 

In summary the procedure is as follows: first check to see if pi and pi 
are both negative at the same point. If so, then the test is failed; if not, 
proceed to determine the intervals where |pi(0) > O,p 2 (0) < 0) and where 
(pi(0) < O,p2(0) > 0}. The test is passed (i.e., no overflow oscillations exist) 
if and only if the maximum of p\(8)/p 2 (0) in the first interval is less than 
the minimum of pi(d)/p 2 (0) in the latter interval. 

In the second method we finely discretize the interval for 8, [0,7r], and 
evaluate pi(0) and p 2 (6) at all the discrete points [0j\. Testing for the 
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hypothesis of Theorem 1 then amounts to testing for the consistency of 
the following set of linear inequalities 



{aiOC2, lo i P2 (e l ) p 2 (e 2 )---P2Mi 



> (1 1 ... 1) (41) 



(There is no loss of generality in assuming that the right-hand side is as 
specified above.) There are standard procedures 20 - 21 for testing for 
consistency of linear inequalities. In any case, Phase 1 of any commer- 
cially available linear programming package does precisely this. If a linear 
programming package is used then the following (dummy) functional 
may be used in the program: minimize (a\ + a 2 ). 

The above method is easily adapted to the generalization of the cri- 
terion which is developed in Sec. VI. 

III. EXAMPLES 

We consider three classes of examples in some detail. In each case we 
tested the criterion by following the second method outlined above. We 
used a linear programming package (written in machine language) made 
available to us by A. M. Odlyzko; the interval [0,7r] was subdivided into 
100 intervals. In every case the computation time was of the order of a 
second. 

3. 1 Example 1: third-order system with repeated real roots 

In this class of examples we take the coefficients to depend on a real 
number p, \p\ < 1, in the following manner: 

ai = -3p, o 2 = -3p 2 , a 3 = -p 3 (42) 

A third-order system with the above coefficients corresponds to an 
underlying linear system with characteristic polynomial (X + p) 3 , i.e., 
the linear system possesses three real roots all repeated at —p. In the 
investigation reported in Ref. 1 we found this class of systems to be in- 
teresting for various reasons. Also, for |p| close to 1 the behavior of the 
system is to some extent representative, at least with respect to oscilla- 
tory behavior, of low-pass systems and high-pass systems, depending 
upon whether p is negative or positive respectively. 
In Ref. 1 we showed for system (1) that 

\p\ > 0.858 => period-3 oscillations exist (43) 

Tsypkin's criterion and the circle criterion (see Sec. 5.2) give 

\p\ ^ 0.500 => no overflow oscillations exist (44) 

An application of the new test yields 

|p | < 0.785 => no overflow oscillations exist 
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Thus, in this class of examples the new test makes a substantial contri- 
bution in reducing the indeterminate region to 0.785 < \p\ < 0.858. 

3.2 Example 2: fourth-order system with repeated real roots 

The class of examples considered here is a natural extension to a higher 
order, m = 4, of the class considered in the previous example. Again all 
the coefficients are determined by one real parameter p where \p\ < 
1: 

fll = -4p, a 2 = -6p 2 , a 3 = -4p 3 , a 4 = - P 4 (45) 

Thus in this example the underlying linear system possesses four real 
roots, all repeated at —p. 

By examining the natural set of four equations associated with a pe- 
riodic solution of period 4, see eq. (12), it is easy to see that a periodic 
solution with elements (1,1,-1,-1) exists if and only if 

(a 4 -a 2 )^l+ |a!-a 3 | (46) 

Thus, we find on substituting for the a's that 

|p| ^ 0.669 =* period-4 oscillations exist (47) 

Tsypkin's criterion and the circle criterion give 

|p| < 0.384 => no overflow oscillations exist (48) 

Application of the new criterion gives 

|p| < 0.610 =» no overflow oscillations exist (49) 

Thus we find that in this example too the new criterion makes an ef- 
fective contribution in determining the region of stability. 

3.3 Example 3: fourth-order filter for sample rate conversion 

The example we consider now, a fourth-order system, was designed 
originally for interpolation and filtering for a terminator in a local digital 
switch. 4 We have reported previously 1 that in its original form the filter 
using saturation arithmetic sustained overflow oscillations. Here we vary 
one of the parameters in the design in order to estimate the modification 
required to guarantee the absence of oscillations. We find that the req- 
uisite variation is large. However, in the process we obtain a measure of 
the effectiveness of the new criterion. 

The example we consider has two pairs of complex poles 

Ai, 2 = Pie ± " h , X3,4 = P2e ±/fl2 (50) 

(The coefficients of the system are not of much interest; however, they 

HIGH-ORDER FILTER USING SATURATION ARITHMETIC 1691 



may be obtained from the information given below.) Also 

Pi = 0.786427817, B l = 37.309784226 degrees (51) 

and 

d 2 = 39.675296075 degrees 

We vary p 2 keeping p lt di, and 2 fixed at the above values; in the original 
design p 2 = 0.952851183. 

In (46) we have given a condition for the existence of limit cycles of 
period 4 with elements (1,1,-1,-1). Translating (46) to the present ex- 
amples gives 

P2 > 0.671 => period-4 oscillations exist (52) 

Tsypkin's criterion and the circle criterion give 

P2 < 0.070 => no oscillations exist (53) 

An application of the new criterion gives 

P2 < 0.665 => no oscillations exist (54) 

This is a rather striking example of the effectiveness of the new criteri- 
on. 

IV. AN APPLICATION 

Here we examine a particular eighth-order system* which has been 
used in an applied research project 22 on a TDM/FDM translator. 23 The 
latter, a system for translating between analog frequency-division and 
digital time-division signals, is an extensively studied application of 
digital filtering. The eighth-order system has been designed to function 
as a low-pass filter with a sampling frequency of 8 kHz and a cutoff fre- 
quency of 2 kHz. Our object here is to demonstrate through an applica- 
tion of the new criterion that it is possible to design the filter as a cascade 
of two fourth-order sections both employing saturation arithmetic such 
that no overflow oscillations are sustained in either section. At least as 
far as overflow oscillations are concerned the margin of safety is adequate 
so that small changes in the coefficients due to quantization of coeffi- 
cients, for example, are not going to cause overflow oscillations to appear. 
It should be emphasized that the result here is not a substitute for a 
design study and the structure suggested may well turn out to be unac- 
ceptable on grounds not related to overflow oscillations. 

The system has four pairs of complex poles; the modulus (p,-) and 



* I am grateful to V. B. Lawrence for bringing this system to my attention. 
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argument (±0,) of each pair is as follows: 

Pi = 0.5115846 , 0i = 32.870 degrees 

p 2 = 0.980274196, 2 = 80.828 degrees 

p 3 = 0.75259969-, 3 = 64.482 degrees 

P4 = 0.892679 , 4 = 75.297 degrees 

We group the first and second pairs of poles together to form one 
fourth-order section and the remaining pairs to form the second 
fourth-order section. The resulting coefficients of the two sections are, 
respectively, 

ai = 1.1718731, a 2 = -1.4912153, a 3 = 0.9075846, 

a 4 = -0.2514954 (56) 

ai = 1.1015710, a 2 = -1.6571120, a 3 = 0.7733805, 

a 4 = -0.45135546 (57) 

Both sections pass the new test. For the first section it may be ascer- 
tained that with 

«i = 6.0819413 and a 2 = 0.07538601 (58) 

the hypothesis of Theorem 1 is satisfied. In fact, the polynomial pi(0) 
is itself positive everywhere except at = 0, where its value is 0. However, 
p 2 (0) > 0. Thus, any positive choice of «i and a 2 chosen suitably small 
will satisfy the hypothesis of Theorem 1. 

For the second section (57), a choice of «i and a 2 for which «iPi(0) + 
«2P2(0) >Oforall0is 

«i = .92348761 and a 2 = .16965636 (59) 

Plots of pi(0), p 2 (0), and aip x (0) + a^ 2 (0) for the second section are 
displayed in Fig. 3. 

It is noteworthy that both sections fail Tsypkin's criterion and the 
circle criterion. 



V. SOME IMPLICATIONS OF THE MAIN RESULT (THEOREM 1) 
5. 1 Overflow oscillations do not exist in second-order systems 

It is well known 12 - 13 - 14 that when the order of the system in eq. (1) is 
two, then overflow oscillations are not sustained. The proofs of this result 
are rather special to second-order systems and to the saturation arith- 
metic. On the other hand, there are the frequency-domain criteria 5 " 11 
for stability which are systems-theoretical results applicable to large 
classes of nonlinearities and systems of arbitrary order. However, we may 
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infer from the results in Ref. 8 that these criteria do not give the result 
that all second-order systems are free from overflow oscillations. 

We show that the criterion in Theorem 1 does give the well-known 
result on second-order systems. Our result is given in Proposition 3. [It 
is assumed that \a 2 \ < 1 and 1 — |ai| - a 2 > 0; these relations are 
equivalent to eq. (3), i.e., the underlying linear system is stable.] 

Proposition 3: Let m = 2 in eq. (1). Also let 

ai = (1 + ai - a 2 ) > and a* ■ (1 - a\ - a 2 ) > (60) 

Then, 

a 1 pi(6) + a 2 p 2 (d)>0 for all 6 (61) 

The proof of this result is in Appendix B. The above in conjunction 
with Theorem 1 shows that oscillations are not sustained in second-order 
systems. 

5.2 Tsypkln's criterion and discrete circle criterion are subsumed by new 
criterion 

The object here is to show that the new criterion subsumes both 
Tsypkin's criterion 5 and the discrete circle criterion 11 when the latter 
criteria are used to determine the nonexistence of oscillations in eq. (1). 
The two closely related frequency-domain criteria are identical when 
applied to the system in eq. (1). 

Tsypkin's criterion 5 is as follows in applications to systems like eq. 
(1) where the nonlinear ity F satisfies 

K min ^F(<r)/o^K max for all a (62) 

If 

(i) £ ajz ~J/[1 - K min £ ojz -J] is finite for all \z\ >1 (63) 

Hi) ~ 3^ Re [E oje-W/a - K min £ aj*-W)] > for all 

*» max *» min 

d in [0,2tt] (64) 
then lim x n = 0; in particular, oscillations do not exist. 

n— ►«■ 

In the case of eq. (1) where F is the saturation nonlinearity, 

X min = and K max = 1 (65) 

so that the effective restriction is (64) which reduces to 

m 

1 - £ aj cos j6 > for all d (66) 

7=1 
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When (65) holds the discrete time version of the circle criterion 11 is 
identical to the above condition. 

In Theorem 1 let c*i = a 2 = l k- Then from the defining relations for 
Pi(0) and p 2 (d) in (24) and (30) respectively, we find that 

m 

aipi(0) + a 2 p 2 (d) = 1 - L aj cosjfl (67) 

Thus, as previously asserted, if either Tsypkin's criterion or the circle 
criterion is satisfied, i.e., (66) is valid, then the hypothesis of Theorem 
1 is also satisfied. 



VI. A GENERALIZATION OF THE MAIN RESULT 

The reader will recall that the main result, Theorem 1, is a direct 
consequence of the rather special passivity properties, stated in Sec. 2.1, 
which are implied by the special features of the saturation nonlinearity 
F. Another key ingredient is that the passivity conditions imply ine- 
qualities on quadratic forms involving circulants. We show here that 
many conditions akin to the ones in Proposition 1 are valid by virtue of 
the properties of the saturation nonlinearity. All or some of these may 
be used to augment the passivity conditions used so far so as to obtain 
improved criteria for the nonexistence of oscillations. 

The following generalized passivity conditions exist* for any I > 1: 

e n (x n -x n -,)^0 n = 1,1 + 1,... (68) 

e n (x n + x n -l)^0 n = l,l + l,... (69) 

where \x n } is any solution of eq. (1) and \e n J is obtained from the solution 
through eq. (6). The proof is similar to that of Proposition 1. Thus in 
Proposition 1 we have used only a very small subset (/ = 1) of all the 
above conditions. 

The interesting fact is that each of the expressions in the above con- 
ditions summed over N, where N is the period of any periodic solution 
of (1), is equivalent to a quadratic form involving a circulant. Thus if X 
= (X h X 2 , ■ ■ • ,XnY are the elements of the periodic solution and 
(Ei,E 2 , . ■ ■ ,En) are the corresponding forcing terms, see eqs. (12) and 
(13), then 

N r ™ m "1 

£ E n (X n - X n -i) = X' I - L ajPJ - P"' - E a;P-' + > X 

n=l L /-I j=\ -I 

(70) 



* The generalized passivity conditions are also valid for negative values of / although we 
do not make any explicit use of this fact. 
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N \ m ml 

E E n {X n + X n -i) = X' I - E ajPJ + P"< + E o/P-""' X 

n=l L 7=1 ; = 1 J 

(71) 

for / = 1,2, .... Hence by transforming X to Z where Z = U*X, U being 
the unitary matrix of eigenvectors ofNXN circulants, we obtain for / 
= 1,2, . . . 



and 



where 



and 



E E n (X n - X n - t ) = E \Z k \* Pl (k2ic/N) (72) 

n=l k=l 



E E n (X n + Xn-i) = E \Zk\ VlWw/N) (73) 

n=l k=l 



m 

pi(0) ± 1 - cos IB - E OylcoB ;0 + cos 0' - 0^} (74) 
y=i 



pj(fl) = 1 + cos W - E a, {cos ;0 - cos 0' - 0*} (75) 

Thus pi(0) and P2(^) defined in Sec. 2.3 correspond to Pi(0) and p'i(0) 
respectively in the present notation. 

Certainly the generalized passivity condition in (68) and (69) imply 
that the expressions in (72) and (73) are nonpositive. We thus arrive at 
the following generalization of Theorem 1: 

Theorem 2: If any convex linear combination of 
Pi(0),Pi(0),P2(0),P2(0)» • • • is positive for all in [0,w], then the system 
in eq. (1) does not have any nontrivial periodic solutions. 

In experiments involving fourth-order systems of practical interest 
we have not found the use of the above generalized criterion to make any 
material difference in delineating the stable systems. In these investi- 
gations we used a linear programming package (Sec. 2.5) to apply the 
test in Theorem 2 with up to six polynomials (the leading six polynomials 
of Theorem 2) being used. However, it is quite possible for substantial 
improvements to exist in other cases. 
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APPENDIX A 
Circulant matrices 

For completeness we collect here some of the well-known properties 
of circulants which are used in the paper. The interested reader may refer 
to Muir 16 and Grenander and Szego 17 for further details and applica- 
tions; Ref. 15 concisely lists some of the main properties. 

We let P denote the primitive N X N circulant: 



P = 



Note that 



and that 





1 — 

10 



1 

— 

— 











(76) 



P N = I 



p/ _ pN-l = p-1 



(77) 



(78) 



A polynomial of arbitrary degree in P is a circulant. AnNxN circulant 
C, 



N-l 

C = L cjPJ 

j = 



has as its eigenvalues 



iV-l 



£ cje-^ 2 *'" k = 1,2, ...,N 
;'=o 



(79) 



(80) 



AWNXN circulants have as eigenvectors Uk, k = 1, . . . ,N, given in eq. 
(16). The matrix U with the eigenvectors as columns is unitary, i.e., 

U*U = / (81) 

APPENDIX B 

Proof of proposition 3 

We prove here the assertion in Proposition 3, namely, for second-order 
systems 

q(0) ± (1 + ai - a 2 )pi(0) + (1 - a x - a 2 )p 2 (6) > for all d (82) 

For second-order systems 

Pi(0) = (1 + ai) - (1 + ai - a 2 ) cos 6 - a 2 cos 26 (83) 
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and, 

p 2 (9) = (1 - a x ) + (1 - ai - a 2 ) cos - a 2 cos 20 (84) 

We find upon substitution that 

q(d) = -4a 2 (l - a 2 ) cos 2 - 4ai(l - a 2 ) cos + 2(1 + a? - of) (85) 

First observe that at the corner points q is positive: 

q(0) = 2(1 - ai - a 2 ) 2 > and g(x) = 2(1 + a t - a 2 ) 2 > (86) 

Through differentiation we find that minima ofq(9) occur in the interior 
of the region [0,ir] if and only if 

|ai|<-2a 2 (87) 

and that at a minimum 6, 

cos0=-ai/2a 2 (88) 

Evaluating q at such a point we obtain 

q 0) = " (1 + Q2) [-2a 2 (l + a 2 ) + (4ai - a?)] (89) 

>0 
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